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Abstract: We investigate modifications of tlie Lifsliitz black liole solutions due to 
the presence of Maxwell charge in higher dimensions for arbitrary z and any topology. 
We find that the behaviour of large black holes is insensitive to the topology of the 
solutions, whereas for small black holes significant differences emerge. We generalize 
a relation previously obtained for neutral Lifshitz black branes, and study more 
generally the thermodynamic relationship between energy, entropy, and chemical 
potential. We also consider the effect of Maxwell charge on the effective potential 
between objects in the dual theory. 
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1. Introduction 

Since the AdS/CFT correspondence was proposed by Maldacena [1] holography has 
been a useful tool in studying strongly coupled field theories. Specifically, holography 
proposes a duality between the gravitational dynamics in an asymptotically AdS 
spacetime and a conformal field theory on the boundary. 

In recent years the idea of holographic duality has been developed beyond high 
energy physics to describe strongly coupled systems in condensed matter physics, 
such as quantum critical systems [4, 5, 6, 7]. Such systems are difficult to study using 
traditional methods in condensed matter based on weakly interacting quasiparticles 
and broken symmetry. Quantum critical points have a spacetime scale invariance 
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which provides a strong kinematic connection to the some versions of AdS/CFT cor- 
respondence. This scahng symmetry is based on anisotropic scahng transformation 
between space and time known as Lifshitz scahng 

t XH, r A"V, X Ax, (1.1) 

where z{> 1) is a dynamical critical exponent represents the degree of anisotropy 
between space and time; manifestly z = 1 exhibits relativistic systems. This scaling 
property (noted previously in other contexts [2]) holographically is represented in 
the following form of the spacetime metric [3]: 

ds' = f (^-r^'dt' + ^ + r'dx"^ , (1.2) 

where the coordinates (t, r, x*) are dimensionless and the only length scale in the 
geometry is £. Metrics asymptotic to (1.2) can be generated as solutions to the 
equations of motion that follow from the action: 

/ = y d-+'xy^ (^R-2A- \f^uF^^^ - \H^,m^ - ^B,B^^ , (1.3) 

where A is the cosmological constant, = d^^A^] with representing the Maxwell 
gauge and if^^, = d[^Bi,^ is the field strength of the gauge field B^ with mass rn? = C. 

Recently a similar action to (1.3) but without a massless gauge field (the Maxwell 
field) has been introduced in four dimensions [8, 9], such that background metrics 
with the anisotropic scale invariance (1.2) are obtained as exact solutions. An exten- 
sion of these results to 5-dimensions including the massless gauge field was carried 
out shortly afterward [10, 12]. The introduction of the Maxwell field introduced a 
new length scale, allowing the gravitational system to undergo phase transitions. 
Working in (4-|-l) dimensions [10], an exact charged solution for z = 6 was obtained, 
and the relationship between black hole temperature and charge was numerically 
computed. Exact solutions in (4-(-l) dimensions [12], and then [n + 1) dimensions 
[13, 14] were subsequently obtained, where z = 2n — 2. These models yield a holo- 
graphic description of a strongly coupled quantum critical point in (n+1) dimensions 
with asymmetric scaling that models the anomalous specific heat found at low tem- 
perature in many heavy fermion compounds [13]. 

While some work has been done to obtain general numerical solutions for charged 
Lifshitz black holes with planar topologies [10, 12, 13], the general character of such 
solutions for general topologies and in arbitrary dimension has not been explicated 
up to now. Here we investigate modifications of the Lifshitz black hole solutions due 
to the presence of Maxwell charge in higher dimensions for arbitrary z. As with their 
neutral counterparts, we find that the behaviour of large black holes is insensitive 
to the topology of the solutions, whereas for small black holes significant differences 
emerge. 



-2- 



We also investigate the thermodynamic relationship between energy, entropy, 
and chemical potential, obtaining a generalization of a relation previously obtained 
for neutral Lifshitz black holes [15, 16, 17]. 

Accordingly, in Sec. 2 we introduce the field equation for the metric functions 
and gauge field while our method is a slightly different from what has been done 
in [13]. In section 3 we discuss the behaviour of solutions at large r and then go 
on to examine near horizon expansions in section 4. Section 5 starts by comparing 
numerical solutions in different dimensions for uncharged black holes and follows by 
investigating the effect that the Maxwell charge has on their behaviour. We consider 
the thermal behaviour of these black holes and how charge affects their temperature 
in section 6. We then go on to consider how charge modifies the potential between 
two particles in the dual theory by investigating corrections to the Wilson loop in 4- 
dimensions. We compute the conserved charge and the relationship between energy, 
entropy, and electromagnetic potential in section 8. We close our paper with some 
concluding remarks. 



2. Field Equations in (n + l)-dimensions 

Using the variational principle the field equations that follow from the action (1.3) 
are: 



du.B 



fj,!/ 



CB 



H, 



(2.1) 
(2.2) 
(2.3) 
(2.4) 



where the equations for the massive gauge field have been rewritten in first-order 
form and where 



1 / I 
2 



1 



-BpE'^g^y — Bi^^By 
(2.5) 



is the energy-momentum tensor of gauge fields. 

The (n + l)-dimensional metric preserving the basic symmetries (1.1) under 
consideration can be written as: 



(2.6) 



3 



where 



dnl 



n—l i—1 
i=2 3=1 

i=l 



k = l 
k = -l 

k = 



(2.7) 



represents the metric of an {n — l)-dimensional hypersurface with constant curva- 
ture (n — l)(n — 2)k and volume Vn-i- The hypersurface is S"'~^, or H"'~^, 
respectively, for = 1, or — 1. 

The gauge fields are assumed to be 

At = ir'K{r), Bt = qir' f{r)j{r), H^, = q£zr'-'g{r)h{r)f{r), (2.8) 

with all other components either vanishing or being given by antisymmetrization. In 
order to get the asymptotic Lifshitz geometry (1.2) we demand /(r) = g{r) = h(r) = 
j{r) = 1 and /t(r) = as r — oo, which in turn imposes the following constraints 

1)^ o 2(z-l) 



C 

A : 



(n 



£2 ' - z 

2P ' 



Applying the ansatz (2.6) to the equation (2.4) yields: 

Q 



-Jg 



(2.9) 



(2.10) 



where Q is an integration constant related to the Maxwell charge (as we will discuss 
later) and we have chosen boundary conditions such that the Maxwell vector potential 
vanishes at the horizon. Substituting (2.9) and (2.10) into eqs. (2.1-2.3), the field 
equations reduce to the system of first order differential equations 



df 



f 



{2[{n-l){z-l)f-ziz-l)h' + {z-iy + niz-2)+n']r'g' 
dr A[n — Ijr"^ 

+2(n - 1) [{n - 2)kfg^ -{n + 2z- 2)r^] - gV2(2-")^2|^ ^2.11) 
t = 4(;7^^2 [in - l)iz - l)f + .(. - l)h^ -iz- If - n[z - 2) - n^] rV 
-2(n - 1) [{n - 2)HV - nr^] + gV2(2-")^2|^ (2.12) 



dr 4(n — l)r2 

+2(n - 1) [{n - 2)kfg 

dh 



^dr ^ ~ ^) 



{2 [{n -l)iz- l)f - z{z - l)h^ + iz- if + n{z - 2) + 

2)r2] - g2^2(2-n)^2| ^ ^^^^ (2.13) 

(2.14) 



[n 



One can easily check that the above system will reduce to the four dimensional 
Einstein case introduced in [8, 9] if one put n = 3 and z = 2 and Q = 0. 
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3. Exact Solutions 



In the previous section we introduced a set of four first order differential equations 
which in general can not be solved analytically because of non-linearity. However, 
it is possible to find some exact solutions under certain assumptions. Indeed, if one 
chooses the fields such that 



fir) = j(r) 



9{r) 



(3.1) 



then the set of four ODE's reduces to three. One may easily solve the differential 
equation for h{r) which yields: 



h{r) 



1 + 



h 



^n—l ■ 



(3.2) 



where h is an integration constant. Now, considering the assumption (3.1) and 
substituting eq. (3.2), one may find that Eqs. (2.12) and (2.13) would be satisfied 
provided: 

2r^(2 -2n + z) + 2z{z - l)h'^] + 4zr"+^(2z2 + 2n - z{n + 1) - 2)h 

-4HV"+^(n - l)(n - 2){z - 2) = 0. (3.3) 

There are two possibilities for solving this equation. 

First, we can set z = 2{n — 1). Then the two first terma in (3.3) are eliminated. 
By solving the remaining terms for h we have 



{n - 2)2HV"-3 
2{n - l)(3n-4)' 



h 



(3.4) 



Since h should be a constant then the solutions are consistent if /c = or n = 2, 3. 
The exact solutions for /c = and arbitrary n are given by [11]: 



h 



f = f 



1 



2(n - l)V2'^-2^ 
where there is just one horizon located at 



K 



Q 



{n - l)r"-i 



1 



(3.5) 



2{n 



l/(n-l) 



(3.6) 



and the metric functions start from zero at the horizon and monotonically increase 
to approach unity at infinity. 

For n = 2 the metric is the same for all values of k and the solution is 
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Figure 1: Metric functions /(r) in 4-dimensions with: Left) Q = 0.5 for k = —1 (green), 
/c = (red) and A; = 1 (blue). Right) Q = 30 for ah /c's. 



while in 4 dimensions {n = 3) they differ for different /c's [14]: 

^20r2' ^ ^+i0r2 400r4 Sr^' " 2r2 r^y* ^'"'''^ 

One may note that in contrast to RN black holes, the four-dimensional metric func- 
tions (3.8) have just one horizon at: 

ke' ^AkH* + 5og2 

Tn = h -^^ . 3.9 

° 20 20 ^ ' 

For small charges - which means small black holes - the distinction between solutions 
with different k but the same charge are manifest. However as charge increases the 
black hole radius becomes larger and the distinction between metrics of different k 
with the same charge becomes very small, as shown in Fig. (1). 

The other alternative occurs for < 2; < 1. In this case the first two terms in 
eq. (3.3) would be canceled if one chooses: 

Q"^ = 2z{l - z)h^ . (3.10) 

Now solving the remaining terms in (3.3) for h yields: 

(n-l)(n-2)(2-^)/c£V-3 

z{2 + z-2z'^)-z{2-z)n ' ^ ' 

which is a constant for charged solutions if n = 3. Substituting Eqs. (3.10) and 
(3.11) and solving the field equations for n = 3, the solutions are given by: 



^(^2-2^ + 2)r2' ' ' {z^-2z + 2)r^' 

^.i3^AU-"S) (3.12) 



2(^2 -2^ + 2)2 r2 \^ r^-2 
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Therefore the 4 dimensional metric function with < z < 1 is a unit function for 
k = 0, a naked singularity for k = 1 and a black hole with one horizon located at 
rl = e/{z^ -2z + 2) for k = -1. 



4. Solutions at large r 

In the previous section we wrote down the field equations to higher dimensions and 
arbitrary z (equations (2.11)-(2.14)), and here we review the general form of the 
asymptotic behaviour of the solutions [13]. 

We begin with linearizing the system in (?7, + l)-dimensions. Since we require the 
general metric (2.6) to asymptotically approach the Lifshitz one (1.2), we investigate 
the behaviour at large r by applying a small perturbation for the fields 

fir) = 1 + wfi{r), 
g{r) = 1 +wgi{r), 
jir) = 1 +wji{r), 

h{r) = l + whi{r). (4.1) 



In the charged case we must also consider constraints on the behaviour of the 
gauge field to ensure Lifshitz asymptotics. We first note that the quantity Q is 
proportional to the electric charge per unit volume, since 



167rfifc Js 167rfife J 167rl]fc J r"^-^ 167r 

(4.2) 

where and are the unit timelike and spacelike normals to a sphere of radius r. 

Since fir) does not contribute to the equations for g{r\ hir) and jir) we can first 
study the set of equations involving {(?, Inserting the perturbative expansion 
(4.1) into equations (2.12-2.14), we obtain the equations for small perturbations 




— n 


z{z-\)l{n- 1) 


z 




n — 1 


\ — n 


n 




n + 2 - 2) 


z{n + z- 2)/{n - 1) 


1 - 






+ ( 91 +(!^)|0| (4.3) 

where we have included the Maxwell gauge field as a first order perturbation which 
means that we substitute Q"^ / r'^^"'~^^ = wQ'^ / r'^^"'~^\ since its falloff may be slower 
than other terms in the metric functions. We have also rescaled r r/i. Note that 
for 7^ a universal 1/r^ mode also contributes to the metric functions. 
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A detailed discussion of the large r expansion is given in the appendix. The 
eigenvalues of the matrix of coefficients is obtained via straightforward calculation 



z + n 



1 



\ 



z+n— 1 
2 



+ 



y/9z^-2{3n+l)z+{n'^+6n-7) 



(4.4) 



2 



v 



z+n—1 
2 



A/922-2(3n+l)z+(n2+6n-7) 



2 



indicating there are three independent eigenmodes. 

The two upper modes in (4.4) are always decaying. However the lowest mode, 
depending on the value of z, will either be a growing mode if 2: > n — 1, a zero mode 
(independent of r) if 2 = n — 1, or a decaying mode if 2; < n — 1. Then if one looks for 
the solutions that asymptotically approach a Lifshitz fixed point, one has to discard 
the lowest mode in (4.4) for z > n — 1 by demanding the amplitude of this mode 
vanishes at large r. 

Removing zero or growing modes can be done numerically in both uncharged 
or charged solutions by fine-tuning initial values upon solving the field equations. 
It is obvious that one has to remove the growing mode for z > n — 1. The zero 
mode for 2 = n — 1 becomes a marginally growing mode when non-linear corrections 
are included [8]. Such a mode with positive amplitude cannot approach the Lifshitz 
geometry (1.2) asymptotically and therefore must be removed. On the other hand, 
a zero mode with negative amplitude will slowly decay; in this case a consideration 
of the evolution of /(r), (ignored thus far), indicates that it goes to zero as r — )■ 00 
again yielding unacceptable asymptotic behaviour. 

For 2; < n — 1 all three eigenmodes decay, and so no fine-tuning of initial values 
is required if one simply wishes to obtain solutions that asymptote to the Lifshitz 
metric (1.2). For a given event horizon size we have a family of solutions that are 
all asymptotic to the Lifshitz metric (1.2), but with different fall-off rates. By an 
appropriate choice of boundary terms, such solutions will have finite energy [20, 21]. 

5. Near horizon expansion 

In order to investigate the near horizon behavior of the solutions, consider the fol- 
lowing expansions for the fields: 



f{r) = foVr - ro(l + /i(r - Tq) + /2(r - rof H ), 

9{r) = 1 ^° (1 ^g\[r - ro) + g2{r - r^f ^ ), 



r — ro 



j(r) = jfi^r - ro(l + ji(r - tq) + nij - r^f H ), 

h{r) = ho{l + hi{r - ro) + /i2(r - tq)^ + ■■■), 
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f(r) 0.6- 
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Figure 2: Metric functions /(r) for zero modes in 4, 5 and 7-dimensions from bottom to 
top respectively. Left) rg = 20 for all k. Right) tq = 0.6 for k = 0. 



where ro is the horizon radius. All other constants in the series solutions (5.1) can 
be obtained in terms of ro and ho as well but we omit them here because they are 
too lengthy. 

Inserting this into eqs. (2.11-2.14) and demanding the coefficients for each power 
of r — ro vanish determines all constants in terms of ro and Hq. For example go is 
given by 



^0 =^^0 



1/2 



2(n - 1) 



rt'''iQl-Q')-2ziz-l)hl 



(5.2) 



where 



Q 



c — ^4-2n 



{ [{z - 1)2 + (2 - 2)n + n^] rl + {n- l){n - 2)kf} . (5.3) 







Requiring qq to be real we obtain the constraint 



\Q\ < JQI - 2z{z - l)hli 



,2n 



(5.4) 



generalizing a similar constraint obtained for Q = black holes [8, 9]. Indeed, for 
Q = the constraint (5.4) reduces to 

[{z-lf + {z-2)n + n^ - z{z — l)h'^ rQ + (n — l)(n — 2)ki'^ > (5.5) 

imposing either an upper bound on |/io| if = 0, 1 or imposing a lower bound of 

/ (n-l)(n-2) -r J . 



For Q ^ equation (5.4) is more usefully written as 



[z - 1)2 + {z- 2)n + - z{z - l)/i, 



0„2 



{n - \){n - 2) 



r^ + kt > 



2£2(n-l)(n-2)r, 



2(n-l) 




(5.6) 
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Figure 3: Metric functions /(r) for zero modes with ro = 0.6 in 4 (red), 5 (blue) and 7 
(green)-dimensions for: Left) k = 1; Right) k = —1. 



or alternatively as aX + k > where 

{z - 1)2 + iz- 2)n + ^2 - z{z - l)hl 



a 



{n - l){n - 2) 



Q2£-2(n+l) 

2(n-l)(n-2) 



(5.7) 



upon setting Tq = X£'^ iy 2{n-i)\n-2) ) " " • ^^^^ analysis of this equation involves 

finding the intersection points of a straight line in X with a curve behaving as X^~". 
The line intersects the origin for k = 0, and has an intercept at ±1 for k = ±1. 
For /lo — the line always intersects the X^"" curve, and so any values of X larger 
than this are admissible. Hence for any given Q there is a lower bound on the 
size of the black hole regardless of the value of k. The lower bound is largest for 
k = —1 and smallest for k = 1. As ho increases the slope of the line decreases, 
and so this lower bound increases. For = 0, — 1 it becomes infinite as a — )■ 0, i.e. 



1 

n-l 



/in -> 



(2-l)^ + (2-2)Tl+n^ 



Indeed for small a we find 



X > 



a 



n-l 



[n 



2){n 



[^2n-2 ^ (4n2-8n + 3)(n 



(z-lY +(z-2)n+n'^ 



the 



and we see that the lower bound on X diverges as a — )■ 0. 

However for k = 1 the situation is quite different. As — )■ z(z~i) 
lower bound on X remains finite. As becomes even larger , a larger lower bound 
on X appears since now a < 0. As increases, a becomes more negative, and the 
lower bound on X continues to decrease. Eventually a limit of X > n^-^ is reached 
at which the line is tangent to the curve, where a = (1 — n)n^^ , or 



z{z - l)hl -{z- 1)2 + (^ - 1)n + 
(n- 1)2(^-2) 



g2^n£-2{n+l) 

1{n - l)(n - 2) 



1 
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For values of ho larger than this, no black hole solutions exist; instead there is a 
naked singularity. 

Alternatively if we fix tq then demanding qq to be real implies Ql > (and so 
\Q\ < \Qc\) as well as 



Positivity of Ql is always satisfied if A; = 0, 1. However for = — 1 it is satisfied 
provided tq respects a lower bound as 



The value Q ^ Qc is the extremal limit of the charged Lifshitz black hole [12, 13], 
defined in eq. (5.3); it corresponds to the value of Q for which the temperature 
vanishes. 

6. Numeric Solutions in (n + 1) -dimensions 

In order to find numeric solutions for the system of ODE's (2.11-2.14) we apply 
the shooting method by adjusting some initial values for the fields f{r), g{r), j{r) 
and h{r). To find the initial values we employ the near horizon expansions (5.1) by 
substituting tq + e for r where e <^ 1 , choosing some values for tq and ho such that 
after applying numerical method the metric and field functions approach unity at 
large r. 

6.1 Uncharged Einstein-Lifshitz Solutions 

In order to study uncharged Einstein-Lifshitz solutions we put Q = in ODE's (2.11- 
2.14) and use the numerical method to find the solutions. We already explained there 
are different situations depending on the value of critical exponent z: 

I) We begin with an examination of situations where the zero mode is present, 
which implies z = n — 1. For a given horizon radius tq the solution is unique, 
requiring a fine-tuning of initial values to approach the Lifshitz metric (1.2) at large 
r. We find that the solutions are similar to the n = 3 case [9]: for each n, when 
k = the same value of Hq yields black hole solutions for all values of tq. Specifically, 
hf) = 1.3737, 1.3344, 1.3084 in 4, 5 and 7 dimensions respectively. For k = 1 (—1) it is 
necessary to systematically adjust the value of Hq upward (downward) relative to the 
value for k = 0. Also, for each n large black holes are almost indistinguishable for 
different choices of k, see Fig. (2-left). Distinctions start to appear between solutions 
of the same dimension but different k as the horizon radius gets smaller. 

We can also probe the behaviour of our solutions as a function of dimensionality. 
Looking at the left-hand side of Fig. (2), we see that for large tq, where the solutions 




(5.8) 




(5.9) 
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r r 

Figure 4: Metric functions /(r) for A; = — 1 in 7-dimensions with z = 2: Left) rg = 
0.92. Right) ro = 20 while the massive gauge field strength /iq increases and temperature 
decreases from top to bottom in both cases. 

are almost the same different values of k, as the dimension increases the solutions 
approach (1.2) more rapidly as r increases. However the overall behaviour of the 
solutions is essentially the same: the metric function /(r) starts from zero at the 
horizon radius and monotonically increases to asymptotically approach unity at large 
r. The situation is the same for small black holes when k = 0, shown at the right in 
Fig. (2). 

However, for small black holes with k = 1 or k = —1 the situation differs, as 
depicted in Fig. (3). For k = 1 in 4-dimensions, the metric function /(r) starts 
from zero at the horizon radius increasing to reach a maximum above unity, then 
decreasing to approach unity at large r. However in dimensions greater than 4, 
the metric function /(r) starts from zero at the horizon radius and monotonically 
increases to asymptotically approach unity as r goes to infinity. For k = —1 the 
situation is completely vice- versa. Examples are shown in Fig. (3), with /c = 1 on 
the left and k = —1 on the right. 

II) z < n — 1: For the series solutions at large r all three eigenmodes resultant 
from the small perturbation are decaying, and so no fine-tuning is required. In this 
case it is possible to find solutions that all approach unity at large r by fixing tq and 
then searching for a value of Hq such that the boundary conditions are satisfied. We 
find that ho is not unique; in other words there is a continuous spectrum of values of 
ho that all respect the boundary conditions. Numerical exploration shows that for 
strong massive gauge field, i.e. a larger value of ho, the metric function /(r) starts 
from zero at the horizon radius and monotonically increases to asymptotically reach 
unity. However, for a smaller value of ho, f{r) starts from zero at the horizon radius 
and as r increases, passes through a maximum value greater than unity, and then 
approaches unity as r goes to infinity. Moreover, black holes with stronger massive 



- 12 - 




1.2 -| 1.2-| 




Figure 5: Metric functions /(r) in 4, 5 and 7-dimensions with z =4, 5 and 7 from bottom 
to top, respectively. Left) rg = 20 and all values of fc's. For a given dimension the different 
values of k lie almost exactly on the same curve. Right) tq = 0.6 and /c = . 




oH — • — I — ■ — I — • — I — ■ — I — • — I — ■ — I — • — I — ■ — I — • — I oH ■ 1 ■ 1 • \ ■ 1 

0123456789 1 2 3 4 



Figure 6: Metric functions /(r) with tq = 0.6 in 4 (red), 5 (blue) and 7 (green)-dimensions 
and z = 4, 5 and 7 respectively. Left) k = \. Right) k = —1. 



gauge field are colder than those with weaker ones. For example, in 7 dimensions 
with z = 2 and k = —1 then, for a fixed tq, there is a family of solutions that have 
qualitatively different behaviour compared to each other but all are asymptotic to 
the Lifshitz background at large r. An illustration of this is in Fig. (4) for both small 
and large black holes. We find the same situation in other dimensions with different 
choices of z and k provided z < n — 1. 

Ill) z > n — 1: In this case one eigenmode amongst three is growing. To have 
solutions asymptotic to (1.2) the initial values of the parameters must be adjusted 
such that the amplitude of this mode is removed otherwise the fields diverge so fast 
at large r. Compared to the zero mode the fine-tuning is numerically more delicate 
- approaching (1.2) at large r needs to more accuracy in digits of ho. However, the 
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Figure 7: Metric functions /(r) for tq = 20, z = 2 and all k. Left) in 5 dimensions 
and (5 = 0, 30000, 46600 from top to bottom [Qc ~ 46600]. Right) in 4 dimensions and 
(5 = 0, 1100, 1780 from top to bottom [Qc ~ 1788]. In each case /c = is red, /c = 1 is 
blue and A: = — 1 is green; these different cases are just barely distinguishable for such large 
values of r . 

results are similar to the zero mode case: for each n large black holes are almost 
indistinguishable for different choices of k and the metric function /(r) reaches unity 
faster as the dimension of spacetime increases, as shown in the left side of Fig. (5). 
Distinctions start to appear between solutions of the same dimension but different k 
as the horizon radius ro gets smaller, illustrate in the right side of Figs. (5) and in 
Fig. (6). 

6.2 Charged Einstein-Lifshitz Solutions 

We now consider how the Maxwell gauge field changes the solutions. First, the 
requirement [13] that the field strength [r^n)' of the massless gauge field vanishes as 
r goes to infinity implies z < {n — 1). For the zero mode case {z = n — 1), we find 
that ii Q <^ Qc these solutions are not significantly different from uncharged ones. 
However distinctions start to appear as the Maxwell charge Q becomes comparable 
to Qc- Indeed in each dimension, for large black holes both Qc and the metric/gauge 
functions are almost indistinguishable for differing values of k and the same charge, 
as shown in Fig. (7). This is not the case for small black holes. Here Qc significantly 
depends on k and the solutions are distinguishable for differing k. However in general, 
as Q increases the solutions either approach unity at large r less rapidly or become 
less sharply peaked for intermediate values of r in situations where such peaks exist. 
For example, three different solutions for small black holes of the same radius but 
different Maxwell charges are illustrated in Fig. (8). The left-hand side of the 
diagram is for A; = 1 in 4-dimensions and the right-hand side is for k = —1 in 
7-dimensions. 
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Figure 8: Metric functions /(r) for different dimensions. Left) 4-dimensions with tq = 0.6 
for zero mode and A; = 1 for Q = 0, 1, 2 from top to bottom, respectively [Qc ~ 2]. Right) 
7-dimensions with rg = 0.6 for zero mode and k = —1 for Q = 0, 0.3, 0.41 from up to 
down, respectively [Qc ~0.417]. Increasing values of Q are coloured red, blue, and green 
respectively in each case. 




Figure 9: metric functions /(r) in 5 dimensions for ro = 0.6 and z = 2 for: Left) k = 
and Q = 0, 0.8, 1.259 [Qc ~ 1.259], Middle) A; = 1 and Q = 0, 1.1, 1.77 [Qc ~ 1.77],Right) 
k = —1 and Q = 0, 0.1, 0.176 [Qc ~ 0.176] from top to bottom, respectively. Increasing 
values of Q are coloured red, blue, and green respectively in each case. 

If z < n — 1, for a fixed tq there is a family of solutions depending on fall-off 
rate. Consequently to investigate how the Maxwell charge modifies the solutions 
by demanding the same fall-off rate (for a given value of k) we search for possible 
charged black holes. Numerical exploration reveals that with a fixed ro, we cannot 
find a solution if Q > Qc, where Qc is given by eq. (5.3). Fig. (9) shows /(r) for 
small black holes in 5 dimensions with z = 2 and k = 0, 1, —1 for different Q. 

Corresponding plots can be constructed for g{r) and h{r)] we have not repro- 
duced them here. 
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Figure 10: Temperature vs rg with Q = for zero modes in 4(red), 5(blue) and 7(green)- 
dimensions for: k = —1 (dash-dots), k = (dots) and k = I (sohds). 

7. Thermal behavior 

The temperature of Lifshitz black holes is easily evaluated using standard Wick 
rotation methods, yielding the result 



where Qq is given by (5.2) and /o is determine so that the metric has asymptotic 
behavior given in equation (1.2). The temperature clearly depends on the spacetime 
dimensionality {n + 1), the critical exponent z, the topological parameter k, the 
horizon radius tq, and the Maxwell charge Q. 

The behavior of temperature versus tq for uncharged solutions Q = and zero 
modes z = n — 1 in different dimensions is depicted in Fig. (10) for k = —1,0 and 
1. For all topologies we see that the temperature increases with dimensionality for 
large black holes, whereas it decreases with dimensionality for small black holes. This 
effect is most pronounced for k = 1, and least so for /c = — 1. In each dimension the 
temperature decreases as the horizon radius shrinks, i.e. small black holes are colder 
than large ones. For k = —1, there is a lower bound on the radius of the black hole. 

For Q ^ 0, since the metric functions /(r) and g{r) are not altered signifi- 
cantly when Q is small compared to its upper bound Qc no significant change in the 
behaviour of the temperature is expected. This upper bound value increases with 
increasing horizon radius. Indeed, with a fixed horizon the temperature goes to zero 
as Q — )■ Qc- We therefore consider the behaviour of temperature as a function of Q 
for fixed horizon radius vq. 

Figure (11) illustrates the situation in 4-dimensions for both small {tq = 0.6) 
and large {tq = 20) black holes. For all topologies the temperature decreases with 
increasing Q as expected, with temperature being largest for the spherical case and 
smallest for the hyperbolic one. For large black holes this latter effect is almost 
completely indistinguishable, whereas for small black holes it is quite pronounced. 



T 




4vr5fo ' 




Figure 11: Temperature vs Q in 4-dimensions with z = 2 for: ro = 0.6 for k = —1 (green), 
= (red) and k = 1 (blue) on the left and ro = 20 on the right; in the latter case the 
curves are almost identical for each k. 




Figure 12: Temperature vs Q in 5-dimensions with z = 3 for: tq = 0.6 for k = — 1 
(dash-dot), A: = (dot) and k = 1 (solid) on the left and ro = 20 for all k on the right; in 
the latter case the curves are almost identical for each k. 

Furthermore, this effect does not hold as dimensionahty increases. In Fig. (12) we 
plot temperature versus Q in 5-dimensions for both small and large black holes. In 
the former case we see that for sufficiently large Q flat topologies yield hotter black 
holes than spherical ones. In the latter case the distinctions between topologies are 
very small - they exhibit nearly the same temperature for different k, a feature 
we find true in all dimensions we have checked. In any dimension the temperature 
approaches zero as the Maxwell charge approaches its extremal value. 

8. Wilson loop 

In 4 dimensions the action for the dual theory is conjectured to be 

S = J dtcfx (^02 - K(VV)') (8.1) 

and if we write V^^ = V x E, where Ej = ejkV^(f), then the boundary action 
(8.1) could be regarded as a gauge theory in (2-1-1) dimensions with a dimensionless 
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Figure 13: Boundary length L vs in 4-dimensions with z = 2 and tq = 0.5 for k = 1 
(blue), A; = (red) and k = —1 (green) while in each case dot is for Q = 0, dash-dot for 
Q = Qc/'i and solid for Q = 2Qc/3. 



coupling constant [8, 7]. We introduce Wilson loops by joining charged particles 
on the boundary that are connected together in the bulk via a string. These loops 
contain information about the force acting between particles charged under the gauge 
fields in the dual theory. 

The Euclidean action of this string for a rectangular Wilson loop is the same for 
all values of k and is given by [7, 19] 



5 = ^ y dtdT^det[gABd^X^d,XB] 



27ra' / V ^ \dd 



taking a = 9 and zr = t in the static gauge, with Euclidean time interval A. 
If one extremizes the action one may obtains a constant of the motion as 

j2^2z+2 



p^2.+2 + pg2^2z-2 (I) 



f{rm)r, 



z+l 
m 



(8.2) 



13) 



where > rg is the location of the midpoint of the string for which ^ \r„,= 0. From 
the above expression one may get the boundary length as 

dr 



dd = 2 



9 



^ , 2 / X 22+2 

friz 



4) 



and the regularized potential energy between the two particles 



V 



S _ i 



-oo 

2 / dr 



\ 
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Figure 14: Potential V between two particles as a function of string midpoint in 4- 
dimensions with z = 2 and = 0.5 for A; = (red), k = 1 (blue) and k = —1 (green) for 
Q = (dot), Q = Qc/^ (dash-dot) and Q = 2Qc/3 (solid). 

where fm = f{rm)- 

We plot in Fig. (13) how the behaviour of the Wilson loop for the neutral case 
[8, 9] is modified for nonzero Q. For a given we find that the boundary length 
increases and the potential between objects in the dual theory decreases (Fig. (14)) 
with increasing Maxwell charge. 



9. Conserved Charge 



To explore the conserved charge we redefine the metric and gauge fields as follows 
while we have chosen i = 1 for simplicity through the rest of the paper: 

At = e^M, Bt = ge-^M, Htr = ge^^. (9.1) 

where new fields F, G, i?, J, H can be expressed in terms of the old fields as: 

F{r) = In(rV), G{r) = ln{g/r), R{r) = Inr, (9.2) 
H{r) = \n{zr'^^ghf), J{r) = ln(r^j/), K{r) = In(r^K). (9.3) 

If one inserts the ansatz (9.1) into the action (1.3), one obtains the one dimen- 
sional lagrangian Cid = Cig + Cim with 



l9 



(n-1) 
1 



n 



2F'R' + {n- 2)R" 



^F-G+{n-l)R 



I2 2K-2J-2G 



(9.4) 



where we have used eq. (2.3) to write H{r) in terms of J{r). Applying the same 
method described in [17] reveals that 

C = 2{F' - R')e^-G+in-l)R _ (^2 j/g2J _ ^,^2Ky-F-G+in-m 



„n+z— 1 



[2iz - l)f + 2rf' - q^jizfj + rff + rjf')] - Q 



r K, 



(9.5) 
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is conserved along the radial coordinate r and 



K{r) 



(9.6) 



where D is an integration constant, which is chosen such that K{r) vanishes at the 
horizon, that is 

D = -Qj J^^dr. (9.7) 

We shall evaluate C both at the horizon radius and infinity. Now using (9.6) and the 
large r expansions (11.1), (11.2), (11.3) and (11.8) conserved charge for z ^ n — \ is 
given by: 

2{z-l){n- z-l){z + n-l) ^ 
(n-l) 

and for 2; = n — 1 is: 



C = ^ ,^ ^Ci - QD. (9.8) 



2{n-l) n-V ^ ' 

Using the near horizon expansions (5.1) we find that 

C = ^r^+' = IQirTS, (9.10) 
do 

while we have considered massless field is zero at the horizon and 5" is the entropy 
of the black brane per unit volume which in terms of the horizon radius is given by: 

S = (9.11) 

We remark that the conserved quantity C for z = 1 with g{r) = l/f{r) reduces to 

C = 2r"+V/'+- (9.12) 

(n-2)r"-2' V ; 

which is proportional to the mass of the asymptotically AdS-Reissner-Nordstrom 
solution and gives the metric function as 

/'W = 1 - + 2(n-l)(7?-2)r2("-i)' ^^'^^^ 

which is the well-known asymptotically AdS Reissner-Nordstrom solution with fiat 
horizon. 



10. Finite Action and the Energy Density for Einstein-Maxwell- 
Lifshitz Solutions 

In order to have a finite action we must add some boundary terms to the original 
bulk action. For this purpose, we use the same method proposed in [20] to make 
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[n + l)-dimensional Einstein-Maxwell-Lifshitz(EML) action finite and well-defined. 
We consider / = /buik + -^bdy5 where /buik is given in eq. (1.3) and /bdy is the boundary 
action which for the case of zero curvature boundary can be written as: 



'bdy 



57r 



dM 



K - 



1 



n 



Stt 



hn^F^i/A ~l~ /dcriv ) 



dM 



(10.1) 

where in the first integral, which is added to cancel the divergences from gravity and 
the Lifshitz field, the boundary dM. is the hypersurface at some constant r, hap is 
the induced metric, K is the trace of the extrinsic curvature, K^i^ = V (anp) of the 
boundary (where the unit vector n'^ is orthogonal to the boundary and outward- 
directed) . 

In the grand canonical ensemble we set 5A^ = on the boundary and the 
variational principle is well-defined only if a; = 0. In the canonical ensemble we 
hold the electric charge (given by eq. (4.2)) fixed, which is equivalent to setting 

—hn^F^j^v) = on the boundary [24]. In this latter instance the variational 
principle is well-defined provided we select uj = \. The last term Jderiv in (10.1) is 
a collection of terms involving derivatives of the boundary fields, which could in- 
volve both the curvature tensor constructed from the boundary metric and covariant 
derivatives of Ba- Since the boundary is fiat and the fields are constants for (1.2), 
this term will not contribute to the on-shell value of the action for the pure Lifshitz 
solution or its first variation around the Lifshitz background and therefore we ignore 
it henceforth. As explained in [20], an arbitrary function f[B"Ba) is added to the 
action which is due to the fact that on the boundary BaB°' = —q^ is constant for 
Lifshitz solutions. 

The variation of the action about a solution of the equations of motion is just 
the boundary term, 

SI = — [ d"xv^{n„«5/i°^ - (1 - 2uj)n^'F^JA'' - n^H^JB" (10.2) 
167r Jqj^ 

+/'(fi,5")(2fi,55" + B^BpSh'^P) - ^/(5,fi")/i„^5r^}, 



where 



Now, if one defines 



(10.3) 



Sap 

S^ - 



s. 



M 



-h 



16% 

IGvr 



llap + ^{-B,B^r'/\B^Bp - B^B^h^fs) 
[n^^H^^ - zq{-B,B^')~"^B^] , 



167r 



(1 - 2uj)n^F^ 



(10.4) 
(10.5) 
(10.6) 
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Then the general variation of the action is^ 

5/ = y rxiSai35h'"^ + S^5B'' + (10.7) 

Accordingly we can define a stress tensor complex [20] consisting of the energy 
density S, energy fiux £i, momentum density Vi and spatial stress tensor Uij, satis- 
fying the conservation equations 

dtS + d,£' = 0, dtVj + diV) = 0, (10.8) 

where 

S = 2S\ - SiBt - SijA, r = 2S\ - SlBt - ^VA, (10-9) 

and 



V, = -2S\ + SiB, + S\,A,, V\ = -2S\ + SiB, + S'j^A. (10.10) 
Using the exponential ansatz (9.1) for the metric and gauge potentials, we obtain 



S 



9 



[zq'il - h)jg - 2{n - 1)(1 - g)] + (1 - 2uj)Qr'K 



:io.iii 



for the energy density of the black brane. 

Employing the large r expansions given in Appendix for the metric functions in 
the case of z 7^ n — 1 and 2; = n — 1 to eq. (10.11), we obtain the rather curious 
result that 



Stt (n-1) ^'^ ^ 32w(n^z^l)' ^ ^ ±, Z ^ i 



s 



:io.i2) 



T(n-l) 



i6^C^i + 32^[(3 - 4c^)(^ - 1) Inr + 1] 



z = n 



giving a finite energy density if 2; < n — 1 regardless of the choice of u. We re- 
mark that for arbitrary z choosing a; = | yields finiteness of energy. However it 
is straightforward to show that S^^ is divergent for large r, rendering the variation 
(10.7) ill-defined, unless u = 1/2. 
Choosing w = | gives 



lOTT Jj^ 



(10.13) 



in 



(TxzV-h K 



dM 



zq 



(n-D- -^^-B^B- + -n'^F^.A^ ] (10.14) 



^As wc shall see, it is only in the canonical ensemble that subsequent equations in this section 
are defined. Note that 5 A" = in the grand canonical ensemble. 
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for the action. Using equations (9.8-9.10) we obtain 

Ci = -, ^r^"""!!' ^ATS + DQ] z<n-l, (10.15) 

[z — l)[n + z — l)[n — z — 1) 

and applying (10.15) we find 

71 — 1 

8= \TS+QD] z<n-l, (10.16) 

n + z — 1 

Since (10.16) is the expression for energy at infinity one may interpret the coefficient 
conjugate to the electric charge as the chemical potential. Hence 

$ = D z<n-l, (10.17) 

which is consistent with the definition of chemical potential, measured at infinity 
with respect to the horizon, given as [22] 

= D z<n-l (10.18) 

where x'^ = 9t is the null generator of the horizon. For z = 1 setting g = 0, 
substituting (9.13) in eq. (10.11), and expanding the result for large r yields the 
energy density of Reisner- Nordstrom black brane [23] : 

£ = ^^TT^m, 10.19 
167r 

where m = C/n is the mass of the AdS black hole. The quantity C is the conserved 
charge given by (9.12), which at the horizon in terms of thermodynamic quantities 
and horizon radius can be expressed as 

Q 



C = 167r 



TS+Q- 



z = l (10.20) 



and the coefficient of electric charge is exactly AdS chemical potential [23] : 



dr Q 



<^ = D = -Q - = ^ z = l (10.21) 

r"-i (n-2)r^-2 ^ ' 

Therefore, the energy density for 2; = 1 is 

71 — 1 

£ = (T5+Q$). (10.22) 

71 

which, by comparing Eqs. (10.16), (10.17) and (10.22), we see generalizes to 

71 — 1 

£= (T5+Q$), (10.23) 

71 + Z — \ 

for asymptotic Lifshitz black branes. 
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11. Concluding Remarks 

Our considerations of Lifshitz gravity coupled to electromagnetism in (n + 1) di- 
mensions have indicated a rich array of numerical solutions that depend on the two 
parameters ho and Q. Solutions exist for all topologies, and an extremal limit, de- 
noted by Qc and defined in eq (5.3) exists. 

The general effect of Q is to cause metric functions to more slowly approach their 
asymptotic values, other parameters being equal. We illustrated this with the metric 
function /(r): it tends to grow more slowly toward its asymptotic value with increas- 
ing charge (as shown in section 6). We also found that increasing charge decreases 
the boundary length of the Wilson loop (for n = 3) while causing the potential V 
between two particles to grow more rapidly with increasing string midpoint length 

We also found that we can extend the thermodynamics of Lifshitz black branes 
[15, 16, 17] to the charged case, provided the parameter z < n — 1. We have ob- 
tained an expression, eq. (10.23), for the energy density in terms of the extensive 
thermodynamic quantities entropy and charge density and their intensive conjugate 
quantities for asymptotic Lifshitz black branes, generalizing the z = 1 AdS case 
(10.22). Extending to values of z larger than this entails a choice of boundary terms 
that renders the variational principle ill-defined. Commensurate results have been 
obtained in other models [25] whose solutions smoothly interpolate between Lifshitz- 
like and AdS-like behaviour as the ratio of T//i (with /i the chemical potential) varies 
from small to large values. 

Further work in this area will involve obtaining a better understanding of how 
the general effects in this paper affect the dual theories associated with asymptoti- 
cally Lifshitz spacetimes, and of how quantum-gravitational corrections can likewise 
modify such effects. 
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Appendix 

Here we are willing to explore more details of the eigen modes at large r which have 
been roughly introduced in eq. (4.4). Here we ignore the universal mode in the set 
(4.3) by setting A; = 0. The complete solution to (4.3) at large r provided z n — 1 
is 

( ,_CiG^__C2G2 CsCh {n-2z)Q^ 

9^V) ^z+n-l ^(z+n-l+^)/2 ^(z+n-l-^)/2 4(^^ _|_ ^ _ 2)(n — ^ — lyr^-'^''' 
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hi{r) = 
ji{r) = ■ 



Co 



{2n-z- 2)Q2 



.z+n-l j.{z+n-l+^)/2 ^{z+n-l-^)/2 4(^^ _|_ ^ _ 2) (7^ — ^ — 1)- 



C2J2 



C3J3 



j.z+n-1 ^(z+n-l+v^)/2 j,{z+n-l-^)/2 4(ri — 2; — l)2^2n-2 ' 



where 



7 = - 2(3n + l)z + {n^ + 6n - 7), 

— 1) _ 2; — 1 _ z — 1 

Gi = 



/' 1^2' '-^2 
z{z-2 + n) 



n — 1 n — 1 

n — 3z + 1 — 1/7 



(11.3) 

(11.4) 
(11.5) 



n — 82; + 1 + ,/7 , 
•^3 = ^77 7V^(ll-6) 



(n-l)2 ' ^ 2(n-l) ' ^ 2(n-l) 

If one inserts the above expression in equation for the small perturbation of /(r), i.e 



r^fiir) = {n-2 + 2z)g,{r)-^^^h^{r) + {z-l)Jl{r)-J^ 
ar n — 1 4(n — ijr"^" ^ 



{n - 2)k 



then again the universal 1/r^ mode decays away at infinity and one obtains 



2r2 ' 
(11.7) 



h{r) 



+ 



C2-F2 



+ 



C3F3 



+ 



j.z+n-1 j.{z+n-i+^)/2 ^{z+n-i-^)/2 4(n - 2; - 1) (n - l)r2"-2 
with the constraint that /i (r) should goes to zero as r goes to infinity and 



^, (11-8) 



z{z — l)(n — 2; — 1) 



z{n + 1) - - 22;2 + 1 + (n - 1)^ 



(n- l)2(n + ^- 1)' 2(n- l)(n-;2- 1) 

z{n + 1) - 77,2 - 22;^ + 1 - (n - 1)^ 



2(n- l)(n-2- 1) 
For 2; = n — 1 solutions arc given by: 

(n-2) (Ci In r + Ca) (3n - 4) Ci 



(11.9) 



^i(^) = 



r.2n-2 



+ 



2(n- 1)2 

(n-2)(lnr)^+^2^-^)^"" ■ " ^ 



hi{r) 



n — 1 



4(n- l)3r2«-2 
+ 



2(n- 1)^ 



]{n — l)r 



2n-2 ' 



(Cilnr + Ca) 

" 2(n - l)r2"-2 4(n - l)2r2"-2 



+ 



(lnr)2 + 



Inr 



+ 



n-1 2{n-iy 



(11.10) 

is) 



. , , (2n - 3) (Ci Inr + C2) (2n - 3) Ci 

JH?^) = 777 ^^^^ h 



Mr) 



2{n - 1)2 r2"-2 
(2n-3)(lnr)2 + 
(3n - 4) Ci 



4(n - 1)3 r2'»-2 



(2n-3)lnr^ 4n - 5 



n-1 



2(n-l)2 



8(n-l)r2"-2' 



[{n- l)lnr + 1] 



(3n - 4)Q2 



4(n - 1)3 r2"-2 ^^'^ ' 8(n- l)3r2'»-2' 

which are relevant to asymptotic Lifshitz background. 



(11.12) 
(11.13) 
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